Central Path Dynamics and a Model of Competition, II

Growth -- the change in number or size -- and adaptation -- the change in quality or structure -- are key attributes of global processes in natural communities, society and economics (see, e.g. Hofbauer and Sigmund, 1988; Freedman, 1991; Young, 1993). In this paper we describe a model with explicit growth-adaptation feedbacks. We treat it in the form of an economic model of competition of two firms (with several departments) on the market. Their size is measured by their capital, and their quality by their productive power (production complexity). It is assumed that the production complexity of a department or firm is a simple function (that is more general than the one considered in Krazhimskii and Stoer, 1999) of its capital. The model works on both the firm level (competition among the departments) and the market level (competition among the firms). The model shows some empirically observable phenomena. Typically, one of the firms will finally cover the market. The winner is not necessarily the firm with the potentially higher maximum productivity. A long-term coexistence of firms may arise in exceptional situations occurring only when the maximum potential productivities (not the actual productivities) are equal. The analysis is also based on the concept of central paths from the interior point optimization theory (see Sonnevend, 1985; and e.g., Ye, 1997)

Approaching the Ground State of a Quantum Spin Glass using a Zero-Temperature Quantum Monte Carlo

Here we discuss the annealing behavior of an infinite-range $\pm J$ Ising spin glass in presence of a transverse field using a zero-temperature quantum Monte Carlo. Within the simulation scheme, we demonstrate that quantum annealing not only helps finding the ground state of a classical spin glass, but can also help simulating the ground state of a quantum spin glass, in particularly, when the transverse field is low, much more efficiently.Comment: 8 pages, 6 fig

Exact Diagonalization Dynamical Mean Field Theory for Multi-Band Materials: Effect of Coulomb correlations on the Fermi surface of Na_0.3CoO_2

Dynamical mean field theory combined with finite-temperature exact diagonalization is shown to be a suitable method to study local Coulomb correlations in realistic multi-band materials. By making use of the sparseness of the impurity Hamiltonian, exact eigenstates can be evaluated for significantly larger clusters than in schemes based on full diagonalization. Since finite-size effects are greatly reduced this approach allows the study of three-band systems down to very low temperatures, for strong local Coulomb interactions and full Hund exchange. It is also shown that exact diagonalization yields smooth subband quasi-particle spectra and self-energies at real frequencies. As a first application the correlation induced charge transfer between t2g bands in Na_0.3CoO_2 is investigated. For both Hund and Ising exchange the small eg' Fermi surface hole pockets are found to be slightly enlarged compared to the non-interacting limit, in agreement with previous Quantum Monte Carlo dynamical mean field calculations for Ising exchange, but in conflict with photoemission data.Comment: 9 pages, 7 figure

Steady-state stabilization due to random delays in maps with self-feedback loops and in globally delayed-coupled maps

We study the stability of the fixed-point solution of an array of mutually coupled logistic maps, focusing on the influence of the delay times, $\tau_{ij}$, of the interaction between the $i$th and $j$th maps. Two of us recently reported [Phys. Rev. Lett. {\bf 94}, 134102 (2005)] that if $\tau_{ij}$ are random enough the array synchronizes in a spatially homogeneous steady state. Here we study this behavior by comparing the dynamics of a map of an array of $N$ delayed-coupled maps with the dynamics of a map with $N$ self-feedback delayed loops. If $N$ is sufficiently large, the dynamics of a map of the array is similar to the dynamics of a map with self-feedback loops with the same delay times. Several delayed loops stabilize the fixed point, when the delays are not the same; however, the distribution of delays plays a key role: if the delays are all odd a periodic orbit (and not the fixed point) is stabilized. We present a linear stability analysis and apply some mathematical theorems that explain the numerical results.Comment: 14 pages, 13 figures, important changes (title changed, discussion, figures, and references added

The dynamical stability of the static real scalar field solutions to the Einstein-Klein-Gordon equations revisited

We re-examine the dynamical stability of the nakedly singular, static, spherical ly symmetric solutions of the Einstein-Klein Gordon system. We correct an earlier proof of the instability of these solutions, and demonstrate that there are solutions to the massive Klein-Gordon system that are perturbatively stable.Comment: 13 pages, uses Elsevier style files. To appear in Phys. Lett.

Quantum Computation of a Complex System : the Kicked Harper Model

The simulation of complex quantum systems on a quantum computer is studied, taking the kicked Harper model as an example. This well-studied system has a rich variety of dynamical behavior depending on parameters, displays interesting phenomena such as fractal spectra, mixed phase space, dynamical localization, anomalous diffusion, or partial delocalization, and can describe electrons in a magnetic field. Three different quantum algorithms are presented and analyzed, enabling to simulate efficiently the evolution operator of this system with different precision using different resources. Depending on the parameters chosen, the system is near-integrable, localized, or partially delocalized. In each case we identify transport or spectral quantities which can be obtained more efficiently on a quantum computer than on a classical one. In most cases, a polynomial gain compared to classical algorithms is obtained, which can be quadratic or less depending on the parameter regime. We also present the effects of static imperfections on the quantities selected, and show that depending on the regime of parameters, very different behaviors are observed. Some quantities can be obtained reliably with moderate levels of imperfection, whereas others are exponentially sensitive to imperfection strength. In particular, the imperfection threshold for delocalization becomes exponentially small in the partially delocalized regime. Our results show that interesting behavior can be observed with as little as 7-8 qubits, and can be reliably measured in presence of moderate levels of internal imperfections

Nonaffine rubber elasticity for stiff polymer networks

We present a theory for the elasticity of cross-linked stiff polymer networks. Stiff polymers, unlike their flexible counterparts, are highly anisotropic elastic objects. Similar to mechanical beams stiff polymers easily deform in bending, while they are much stiffer with respect to tensile forces (``stretching''). Unlike in previous approaches, where network elasticity is derived from the stretching mode, our theory properly accounts for the soft bending response. A self-consistent effective medium approach is used to calculate the macroscopic elastic moduli starting from a microscopic characterization of the deformation field in terms of ``floppy modes'' -- low-energy bending excitations that retain a high degree of non-affinity. The length-scale characterizing the emergent non-affinity is given by the ``fiber length'' $l_f$, defined as the scale over which the polymers remain straight. The calculated scaling properties for the shear modulus are in excellent agreement with the results of recent simulations obtained in two-dimensional model networks. Furthermore, our theory can be applied to rationalize bulk rheological data in reconstituted actin networks.Comment: 12 pages, 10 figures, revised Section II

Generalization of the interaction between the Haar approximation and polynomial operators to higher order methods

International audienceIn applications it is useful to compute the local average of a function f(u) of an input u from empirical statistics on u. A very simple relation exists when the local averages are given by a Haar approximation. The question is to know if it holds for higher order approximation methods. To do so, it is necessary to use approximate product operators defined over linear approximation spaces. These products are characterized by a Strang and Fix like condition. An explicit construction of these product operators is exhibited for piecewise polynomial functions, using Hermite interpolation. The averaging relation which holds for the Haar approximation is then recovered when the product is defined by a two point Hermite interpolation

Monomer-dimer model in two-dimensional rectangular lattices with fixed dimer density

The classical monomer-dimer model in two-dimensional lattices has been shown to belong to the \emph{``#P-complete''} class, which indicates the problem is computationally ``intractable''. We use exact computational method to investigate the number of ways to arrange dimers on $m \times n$ two-dimensional rectangular lattice strips with fixed dimer density $\rho$. For any dimer density $0 < \rho < 1$, we find a logarithmic correction term in the finite-size correction of the free energy per lattice site. The coefficient of the logarithmic correction term is exactly -1/2. This logarithmic correction term is explained by the newly developed asymptotic theory of Pemantle and Wilson. The sequence of the free energy of lattice strips with cylinder boundary condition converges so fast that very accurate free energy $f_2(\rho)$ for large lattices can be obtained. For example, for a half-filled lattice, $f_2(1/2) = 0.633195588930$, while $f_2(1/4) = 0.4413453753046$ and $f_2(3/4) = 0.64039026$. For $\rho < 0.65$, $f_2(\rho)$ is accurate at least to 10 decimal digits. The function $f_2(\rho)$ reaches the maximum value $f_2(\rho^*) = 0.662798972834$ at $\rho^* = 0.6381231$, with 11 correct digits. This is also the \md constant for two-dimensional rectangular lattices. The asymptotic expressions of free energy near close packing are investigated for finite and infinite lattice widths. For lattices with finite width, dependence on the parity of the lattice width is found. For infinite lattices, the data support the functional form obtained previously through series expansions.Comment: 15 pages, 5 figures, 5 table

Study of multi black hole and ring singularity apparent horizons

We study critical black hole separations for the formation of a common apparent horizon in systems of $N$ - black holes in a time symmetric configuration. We study in detail the aligned equal mass cases for $N=2,3,4,5$, and relate them to the unequal mass binary black hole case. We then study the apparent horizon of the time symmetric initial geometry of a ring singularity of different radii. The apparent horizon is used as indicative of the location of the event horizon in an effort to predict a critical ring radius that would generate an event horizon of toroidal topology. We found that a good estimate for this ring critical radius is $20/(3\pi) M$. We briefly discuss the connection of this two cases through a discrete black hole 'necklace' configuration.Comment: 31 pages, 21 figure